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Constraining Sommerfeld Enhanced Annihilation Cross-sections of Dark Matter via 
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In a large class of models we show that the light scalar field responsible for the Sommerfeld 
enhancement in the annihilation of dark matter leads to observable direct detection rates, due to 
its mixing with the standard model Higgs. As a result the large annihilation cross-section of dark 
matter at present epoch, required to explain the observed cosmic ray anomalies, can be strongly 
constrained by direct searches. In particular Sommerfeld boost factors of order of a few hundred 
are already out of the CDMS-II upper bound at 90% confidence level for reasonable values of the 
model parameters. 
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Introduction - Strong evidences support the exis- 
tence of Dark Matter (DM) in the present Universe [1], 
although its actual nature is still missing. Identifying the 
DM is a major challenge for particle physics and has lead 
to a vast literature on extensions of the Standard Model 
(SM), in which many new particles comply with the re- 
quirements that a DM should fulfill. Over the last years, 
many efforts have been dedicated in building models of 
DM which leave signatures via direct and/or indirect de- 
tection. 

A huge excitement in the indirect detection of DM 
took place after the PAMELA collaboration [2] reported 
an unexpected rise in the positron fraction at energies 
from 10 GeV up to 100 GeV. Moreover, the H.E.S.S. [3] 
and Fermi Large Area Telescope [4 (FermiLAT) collab- 
orations reported an excess in the electron plus positron 
flux at energies above 100 GeV up to 1 TeV. If the DM 
is indeed the source of these observed anomalies in cos- 
mic rays, then for a stable DM the current annihilation 
cross-section should be boosted by a factor of O(10 3 ) 
with respect to the freeze-out annihilation cross-section: 
(cdmM) ~ 3 x 10 _26 cm 3 s _1 . An attractive way of get- 
ting a large enhancement of this cross-section without 
affecting the DM relic abundance is to invoke the Som- 
merfeld effect |5|. 

Recently CDMS-II (H) has reported two events in the 
signal region at 1.64 a confidence level (C.L.). The col- 
laboration has conservatively set an upper bound on the 
DM spin-independent cross-section on nucleon. 

Typically indirect and direct detections of DM probe 
different sectors of parameter space of a model and have 
been studied as separate issues in the literature. In 
this letter we consider the broad class of model where a 
light scalar <fi is responsible for the Sommerfeld enhance- 
ment [7Mll|. In these models we show that there exists a 
tight connection between direct and indirect detection of 
DM. Indeed the coupling responsible for the Sommerfeld 
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FIG. 1: A schematic presentation of the role of the light scalar 
field <j> in indirect (left) and direct (right) detection of DM. 



enhancement also appears in the spin-independent cross- 
section on nucleon, in a unique combination with the 
Higgs portal coupling (<f> — H mixing, H being the Stan- 
dard Model Higgs) as shown in Fig. Q, Refs. [TlTiTl"] . 
We demonstrate that from the direct detection bounds 
it is possible to set limits on the Sommerfeld enhance- 
ment that can arise at present epoch. Alternatively, for 
a given Sommerfeld enhancement, the Higgs portal cou- 
pling is strongly constrained [121 fT5] and can be probed 
at next generation direct detection experiments. 

Light Scalar - Higgs mixing - Considering that 
the DM candidate % is a SM singlet fermion with mass 
M x , stabilized by a Z2 symmetry, and that the field <fi is 
a real singlet scalar, the relevant terms in the Lagrangian 
are: 

-CD \ x Yx<f> + ^^H. (1) 

Once H develops a non-zero vacuum expectation value 
(vev) v, H and (f> mix through the trilinear term cx fJ-<pv, 
v = 246 GeV. Due to H — <fi mixing, the scalar field <fi 
is unstable and ultimately decays to SM fields. Since 4> 
is produced in the current DM annihilation, we demand 
m 4> ^ 1 GeV in order to avoid the antiproton problem. 
Moreover, in order not to spoil the Big Bang Nucleosyn- 
thesis (BBN) predictions the thermally generated 4> par- 
ticles should disappear before the onset of the BBN, thus 
requiring > 10 MeV [161 [17] , It is true that for a 
DM candidate with mass in the 100 GeV — 1 TeV range, 
such a low mass scale of (f> may appear somewhat unnat- 
ural, being up-lifted by radiative corrections. A super- 
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symmetric realization could however naturally solve this 
problem. In this letter, without addressing this issue, we 
assume to be stabilized at the C(l)( GeV) scale. 

Considering the Lagrangian in Eq. 0, the connec- 
tion between the direct detection and the Sommerfeld 
enhancement is schematically presented in Fig. (JT|), with 
a key role played by <f>. The trilinear coupling X x X°x4> 
gives rise to an attractive Yukawa potential between 
X particles, thus enhancing the current annihilation of 
XX C ~> 4>4>- The same coupling is also responsible for 
a spin independent interaction of x with the nucleon n 
through the Higgs portal coupling We express this 
coupling in terms of the mixing angle between H and 
4>, 0H<t> ~ M<^ w /( to h); where ran is the physical Higgs 
mass in the SM. The mixing angle is lower bounded 
9h<p <; 1CT 7 / \fm$] GeV by demanding that the life- 
time of (f> should be less then tbbn ~ 1 s [TBI IT7] . The 
scalar field <j) is indeed thermally produced in the early 
Universe and should decay before the onset of the pri- 
mordial nucleosynthesis to avoid dominating the energy 
density of the Universe [12 . For the mass range of <f> 
we consider, the mixing between <f> and H is also upper 
bounded, 9 HlP < 1CT 2 [JS]. 

The terms in Eq. (|T|) are ubiquitous in hidden sector 
models. The same Lagrangian arises in the case of 4> 
being a complex singlet and developing a vev. The vec- 
tor DM and the scalar DM cases are very similar to the 
fermionic example, although the presence of extra cou- 
plings may modify the phenomenology of the model. If 
a dark sector is gauged under a hidden symmetry, ei- 
ther abelian or non-abelian, then the DM may be con- 
stituted by the additional vector gauge boson of the the- 
ory |13[ [19] . Typically, SM fields are singlets under the 
hidden symmetry. For both abelian and non-abelian 
cases, the extra gauge bosons will couple to the light 
scalar <fi through the kinetic term 1(9^ — igHT a A^)(f)\ 2 , 
<t> being non-singlet under the hidden gauge group. The 
fermionic DM case then just translates to the gauge bo- 
son one by replacing X x with the hidden sector gauge 
coupling gn- In the non-abelian case a kinetic mixing 
term between the SM hypercharge and the hidden gauge 
sector can arise from higher order operators. In this case 
interesting signatures arise |12l 120]. In the case the DM 
is a complex scalar field S in the hidden sector, the gener- 
ality of our argument is somewhat weakened due to the 
presence of the additional coupling fnsS^SH^H. The 
DM directly communicates to the SM sector through its 
Higgs portal coupling. Thus there is an additional chan- 
nel in the DM to nucleon scattering, mediated by the SM 
Higgs through i-channel. Due to the presence of fas, the 
connection is more involved with respect to fermionic and 
vectorial cases and will be discussed elsewhere |21| . 

Sommerfeld enhancement - In the present case, 
the Sommerfeld enhancement is provided by the light 
scalar field which acts as a long range attractive force 
carrier between the x particles. For a review on the Som- 
merfeld effect due to a light scalar field, see Ref. [7J, 
while for non-abelian massive vector fields one can see 
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FIG. 2: Iso-contours of the non-perturbative Sommerfeld en- 
hancement, (S e ), as a function of vo/(a x c) and e^. 



Ref. [21] • We define the coupling constant between the 
DM and <j> as a x = X x /(Att). When the Compton 
wavelength ©(m^ 1 ) associated with (j) becomes larger 
than (a x M x ) _1 , the asymptotic plane wave ip associ- 
ated with x gets distorted. The distortion can be com- 
puted by solving the Schrodinger equation for the at- 
tractive Yukawa potential V(r) = — (a x /r) e~ m * r and is 
defined as S e = |^(oo)/V>(0)| 2 [7_1|S2]. This is equivalent, 
in terms of Feynman diagrams, to the resummation of 
the multi-loop scalar ladder contributions, as shown in 
Fig. 0, the left diagram. The boost factor S e is a func- 
tion of the dimensionless parameters = m^/(M x a x ) 
and e v = f3/a x only, with /? = v le \/c, the relative velocity 
between the DM particles. 

By considering S e in a halo, we need to integrate over 
the relative velocity distribution of the DM halo [7J [S3J . 
Assuming an isothermal Maxwellian distribution with a 
mean velocity vq, the averaged value of S e reads: 

(S e ) = -±(^Y r d e v ele^y/^S e (e v ^). 
V 7rV v o > Jo 

(2) 

It is a function of and Vo/(ca x ) only. We neglect 
the truncation of the Maxwellian distribution, since the 
cut-off on the escape velocity does not significantly affect 
the results. Indeed the enhancement drops rapidly with 
increasing velocities. The corresponding iso-contours of 
(S e ) are shown in Fig. For small values of vo/(a x c) 
the boost factor can rise up to 10 4 , while for large 
vo/(a x c) the boost decreases down to 1. 

Galactic cosmic rays, being measured by indirect DM 
search experiments, could arise from the current annihila- 
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tion of DM . As a result PAMELA and Fermi experiments 
give upper bounds on the present DM total annihilation 
cross-section from the measurement of e + /~ and p — p 
fluxes. For a DM mass ranging from 100 GeV to 1 TeV, 
the boost factor can be allowed up to a factor 0(1000). 
Interestingly, as we see below, these boost factors suf- 
fer stringent constraints from direct detection exclusion 
limits. 

Direct detection - The light scalar <j> is the only field 
responsible for the DM scattering on nucleon due to its 
mixing with the SM Higgs. The elastic cross-section on 
nucleon, mediated by (j) through i-channel, is given by: 



Ho 



Attv 2 



(3) 



where ji n is the reduced nucleon-DM mass, /„ = 0.3, the 
effective Higgs nucleon interaction. The behavior of a^ 1 
is driven by For a mass of <f> in the MeV-GeV 

range, this cross-section can exceed the current upper 
bound on the elastic cross-section given by CDMS-II. The 
lighter the <fi is, the smaller is the mixing angle required 
to be compatible with direct DM searches. 

Actually what is measured in a terrestrial DM detector 
is the differential rate of nuclear recoils, integrated over 
the energy range of the experiment. This quantity is a 
function of the inverse mean velocity of the DM particles 
that can deposit a given recoil energy E r . For details 
about the total rate the reader is referred to [53] and the 
references therein. With respect to the parameters used 
to define (S e ), the differential rate can be rewritten as: 



dR tiftm* \\c 0* m 1 



dE r 



Attv y 



F(x,y,z) 



(4) 



The function F depends on the minimum velocity to pro- 
duce a given recoil energy x = v m i n /vo, the observed ve- 
locity y = Uobs/^o and the escape velocity z = v esc /vo. 
The observed velocity takes into account the motion of 
the Earth around the Sun. Typically, for an isothermal 
Maxwellian halo the range of values are 170km/s < vq < 
270km/s [55J and 498km/s < v csc < 608km/s [26\. The 
total rate is sensitive to the values of the escape veloc- 
ity, the mean velocity and the minimum velocity, which 
depends on the DM mass, E r and the nucleus mass. 

We consider only the CDMS-II exclusion limit since 
it is the most sensitive one in the DM range we are in- 
terested in (10 GeV-1 TeV). The XenonlO [57] is most 
sensitive in the range 7-20 GeV. The experimental upper 
bound is obtained with the maximum gap method |28| 
and is given at 90% C.L.. 

Results and discussions - We begin by presenting 
our results for a mean velocity Vq of the DM particles 
in the Earth's neighborhood of 220 km/s and an escape 
velocity of 600 km/s. The (S e ) is therefore a function of 
the coupling X x and m$ /M x . In each graph two suitable 
values of the mixing angle 9u<j> are taken into account, 
depending on the <p mass. 
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FIG. 3: Iso-contour of the Sommerfeld enhancement as a func- 
tion of the coupling A x and the ratio of the masses m < f,/M x . 
The thick black solid and dashed curves are the exclusion 
limits from CDMS-II for a mixing angle 6h4, = 1(T 6 and 
0H4, = 10 respectively, with = 0.1 GeV. In analogy the 
thick red long dashed and dotted lines are for m^> = 1 GeV 
and 9h<p = 10 -4 and 8h</> = 10 -3 respectively. The right- 
hand side of each curve is excluded at 90% C.L. The gray 
region is excluded by BBN, Eq. [6] for = 0.25 GeV. 



In Fig. ||) we draw the iso-contours of the Sommerfeld 
boost factor and the exclusion limits from CDMS-II for 
Tjiffi = 0.1 GeV (black curves) and 1 GeV (red curves), 
as labeled in the caption. For 0.1 < A x < 1.2, boost 
factors ranging from 1 up to 1000 can be obtained, the 
right range to account for the PAMELA and Fermi e +// ~ 
and p — p measurements. However these values of X x are 
strongly constrained by direct detection limits on DM- 
nucleon cross-section. For example, if = 0.1 GeV, 
M x = 100 GeV and Qh^ — 10 -6 then only small values 
of X x are allowed, namely A x < 0.3. Conversely, the max- 
imum mixing angle allowed from direct detection can be 
inferred, down to 10" 5 (10" 3 ) for m, p = 0.1 GeV(l GeV). 

In Fig. Q the iso-contours of the Sommerfeld enhance- 
ment are depicted and overlapped with the bounds from 
CDMS-II, for a fixed DM mass of 100 GeV. We can 
infer the allowed boost factor as follows. We see that 
Qh4> = 1CT 6 still allows (S e ) > 500 for value of > 0.5 
GeV. A mixing angle of 10 -4 excludes the possibility of 
having boost factor larger than 10. Had we increased 
the DM mass towards 1 TeV the direct searches would 
have been less sensitive. This is due to the huge dif- 
ference in mass between the DM particle and the nu- 
cleus. For example, the CDMS-II detector is made of 
Ge (mGc ~ 73 GeV), hence the most sensitive region is 
around 50-100 GeV in mass. 
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FIG. 4: Iso-contours of the Sommerfeld enhancement as a 
function of the coupling A x and the ratio of the masses 
m < i > /M x for a fixed DM mass M x = 100 GeV. The thick 
solid and dashed black curves are the exclusion limits from 
CDMS-II for a mixing angle 8 H4> = 10~ 6 and e H <t, = 10~ 4 
respectively. The region below the curves is excluded at 90% 
C.L. 



The effect of the velocity on the exclusion limits and 
on the Sommerfeld enhancement is better understood in 
terms of the rescaled coupling \ x \Jc/vq. As shown in 
Fig. ([2]), (S e ) is indeed only dependent on the rescaled 
coupling. Variations of vq and u csc only affect the di- 
rect detection rate through the F(x,y,z)/y factor, see 
Eq. Q. This modification is plotted in Fig. |5| for the 
extremal parameters vq = 170 km/s and w csc = 500 km/s 
(red dotted and black long dashed lines) with respect to 
the reference values v a = 220 km/s and t> osc = 600 km/s 
(red dashed and black solid curves). This time we fixed 
TO^ = 0.1 GeV. The change in direct detection sensitiv- 
ity may be interpreted as a different value of the mixing 
angle saturating the current experimental upper bounds. 
Larger values of the mixing angle are indeed permitted. 
For M x = 10 GeV and 9h 4, — 10~ 5 , the whole range of 
X x (red dotted curve) is allowed for the extremal velocity 
parameters, while being upper bounded in the standard 
case (red dashed line). For 8H4, = 10 -6 (black curves) 
and M x = 1000 GeV, A x = 1.2 is still permitted for the 
extremal velocity values, implying viable boost factors up 
to 500 (long dashed), while from the reference case (solid) 
the allowed boost is at most of 200. In other words for 
a given Qh^ the maximum allowed X x and boost factor 
can be inferred. Although the constraints from direct de- 
tection are slightly reduced for smaller vq and v csc , the 
main results still hold. 

Notice that we have implicitly assumed an isothermal 
Maxwellian velocity distribution for describing both the 
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FIG. 5: Iso-contours of (S e ) as a function of the coupling 
A x /^/«o/220 km/s and the ratio of the masses m^/M x for a 
fixed <f> mass of 0.1 GeV. The black solid and long dashed lines 
are for Qu<t> = 10 -6 and two different velocity parameters: 
vq = 220 km/s, v CBC = 600 km/s and vo — 170 km/s, v CBC = 
500 km/s respectively. The same for the dashed and dotted 
red curves with 9h 4, = 10 _ . The right-hand side of the curves 
is excluded at 90% C.L.. 



Earth's neighborhood (direct detection) and the galactic 
DM halo ((S e )). Our results however hold even in the 
presence of clumpy structures in the DM density pro- 
file. In that case, from Fig. ([5]), we see that the Sommer- 
feld boost factor is shifted by an amount \/v'q°/ J -UQ lump , 

where t>Q lump ~ 14 km/s. As a result higher boost factors 
become compatible with the direct detection. 

When the DM interacts with the nucleus at a typical 
mean velocity vq/c ~ 10~ 3 , a multi-loop scalar diagrams, 
as in Fig. ([IJ, is also present. Potentially it can boost the 
direct detection rate as it does for the annihilation cross- 
section. With respect to the case of DM annihilation the 
bound state forms between the x particle and the nucleus, 
when they have comparable masses [55]. The attractive 
Yukawa potential of the light scalar field is proportional 
to A x #ff^/ n m n /(47rt;), due to the mixing with the Higgs 
boson. The effective coupling giving rise to the enhance- 
ment is suppressed with respect to a x by the small Higgs 
nucleon coupling and the mixing angle is a rjrj ~ 10 -8 for 
0h4> ~ 10~ 3 . In term of (S e ), from Fig. this corre- 
spond to large values of vo/(a x c). Thus no sizable boost 
factors are expected in this case. 

We conclude this section with a remark on the relic 
density of the DM candidate x an d its link with the Som- 
merfeld boost factor [30 by assuming that the only inter- 
actions are given in Eq. Then the total annihilation 
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cross-section of \ can be estimated as: 



Ml 



1 



2 J 



(5) 



Since the mixing angle Qn<i> oc fi<p, the second and third 
terms, which give the annihilation cross-section of x par- 
ticles into H(f) and H, are suppressed for small mixing 
angles. As a result the first term dominates, which gives 
the annihilation cross-section of XX ~^ 4"^ through the 
X exchange in t— channel. Hence from the relic abun- 
dance of DM, n BM h 2 ~ 0.1, one infers A x ~ 0.6 for 
M x ~ 1 TeV, which is in the expected range of values 
that give large Sommerfcld enhancement as well as de- 
tectable elastic cross-section on nucleon. However, this 
conclusion does not hold if the x particles dominantly 
annihilate via some other channels in the early universe. 

Even after the freeze-out of the DM, a small amount 
of x pairs continue to annihilate with enhanced rate be- 
cause of the Sommerfeld effect. The increased cross- 
section may disrupt the 4 He and D abundances, which 
are benchmark predictions of the BBN. In our case the 
dominant diagram is the t-channel <f> production, first 
term in Eq. [5j XX ~^ < P ( t > - The tf> subsequently decays 
mainly into muons, due to its mass range. Assuming 
that BR(0 — > is 100%, the photodissociation of 

the Helium and Deuterium abundances leads to an upper 
bound on the Sommerfeld coupling A x [TtJl, [T7] : 



A x < 0.05 x 



M x 
GeV 



where E V - 1S is the energy transferred to the visible sec- 
tor. Notice that this constraint neither depends on the 
mixing angle nor on the <j> mass. For a DM mass of 25 
GeV, the maximum allowed A x is 0.5, while for M x = 100 
GeV A x = 2.5. The BBN is capable to set a lower bound 
on the Higgs portal coupling and moreover for lower DM 
masses to strongly constraint the Sommerfeld enhance- 
ment. In Fig. [3] the gray region denotes the parameter 



space excluded by BBN, assuming = 0.25 GeV. For 
larger DM masses, starting from 50 GeV, the CDMS- 
II experiment sets the stringent bounds on the allowed 
boost factors, while BBN constraints are relaxed in the 
range of values of A x we consider. 

At present epoch DM is annihilating in the galactic 
halo producing charged leptons, with an enhanced rate 
with respect to the freeze-out rate, due to the Sommerfeld 
mechanism. An increased neutrino flux can be detected 
at neutrino telescopes, which will constraint the boost 
factor [3T] together with direct searches. 

Conclusions - In this letter we analyzed the con- 
nection between direct and indirect DM searches in a 
class of model where a light scalar field is added to the 
SM through a Higgs portal. This scalar acts as a long 
range force carrier and yields the Sommerfeld enhance- 
ment for the current annihilation of DM, as required by 
observed cosmic ray anomalies. The crucial point is that 
through cf> — H mixing the light scalar allows the DM to 
scatter on nucleon, at rates possibly exceeding the ex- 
clusion limit of CDMS-II. In such scenarios the direct 
and indirect detections are fully determined only by A x , 
Qu<t> and m<p/M x . We found that a large part of the 
parameter space is strongly constrained in order to rec- 
oncile the current bounds from PAMELA and CDMS-II 
at 90% C.L.. For example using the tf> — H mixing of 
order 10" 6 , = 0.1 GeV and M x = 1 TeV, the CDMS- 
II upper bound does not allow neither a boost factor of 
more than 200 nor A x > 0.75. While the interplay be- 
tween direct and indirect detection is straightforward in 
the fermionic example we focused on, we stress that for 
all models where a light scalar is present and the DM can- 
didate has zero SM hypercharge such an interplay does 
occur. 
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